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Measurement of a s and the /3 function with the DELPHI detector at LEP 

O. Passon a 

a Fachbereich Physik, Bergische Universitat Wuppertal, Gaufistrafie 20, 42097 Wuppertal, Germany 

Event shape distributions in e + e~ annihilation are determined from the DELPHI data taken between 183 and 
207 GeV. From these the strong coupling a s is extracted with several techniques. Together with the results 
from other LEP2 energies and at about mz this allows both, a combined measurement of a s and a test of the 
scale dependence of the strong interaction. Alternatively the renormalisation group invariant (RGI) perturbation 
theory is applied to measure the j3 function of strong interaction. The results are good agreement with the QCD 
expectation and allow to exclude the existence of light gluinos with a mass below 30 GeV in a model independent 
way. 



1. Introduction 

Neglecting quark masses, a s is the only free pa- 
rameter of QCD, the theory of strong interaction. 
Thus its measurement is a task of fundamental 
importance. Additionally, the scale dependence 
of the strong coupling, as governed by the f3 func- 
tion, is among the important predictions of QCD. 

This note presents both, a measurements of a s 
at different LEP2 energies, and a test of the scale 
dependence of the strong interaction. In both 
cases several techniques are applied. 

From the event shapes Thrust, C parameter, 
heavy jet mass, wide and total jet broadening, 
a s is extracted with four different methods: The 
differential distributions at each energy are com- 
pared to predictions in 0(a 2 ), pure NLLA and 
matched ©(ev^+NLLA, folded with fragmenta- 
tion models. Furthermore, a s is extracted from 
the mean values using an analytical power correc- 
tion ansatz. These measurements, together with 
the results from other LEP2 energies and about 
mz, are finally combined into one a s measure- 
ment for each method. 

The (3 function is extracted with two differ- 
ent techniques: From the a s values directly, and 
from the energy evolution of mean values of event 
shapes. The last method is based on the renor- 
malisation group invariant (RGI) perturbation 
theory [1]. 

These analyses and the theoretical background 
are described in more detail in [2,3] and the ref- 



erences therein. 

2. Data and background 

The analysis is based on data taken with the 
DELPHI detector in the years from 1997 to 
2000 at centre of mass energies between 183 and 
207 GeV. Details on the DELPHI detector and 
the event selection can be found in [2] and the 
references therein. In total a number of «10.000 
hadronic events have been analyses. This is less 
than 1% of the statistics gathered at y/s — mz- 
But not only the cross section at LEP2 is much 
smaller than at LEP1. Additionally one has to 
deal with significant background processes. Most 
important are radiative (ISR) events and the pro- 
duction of W pairs (four fermion events). For de- 
tails we refer again to [2]. 



Determination of at s 
distributions 



from event shape 



From the differential distributions of 1-T, C pa- 
rameter, Ml/E% [s , B max and B sum , a s is deter- 
mined by fitting an a s dependent QCD predic- 
tion folded with a hadronisation correction to the 
data. As QCD predictions 0(a 2 ), pure NLLA, 
and the combined 0(a 2 )+NLLA in logi?-scheme 
are employed. 

In [4] it has been shown that fixing the renor- 
malisation scale to [i 2 = E 2 m results in a marginal 
description of the data. Therefore, the experi- 
mentally optimised scales from [4] are used for 
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the 0{a 2 s ) fits. For the NLLA and the combined 
NLLA+0(a^) fits, \i is set equal to E cm . 

3.1. Theoretical uncertainties 

All perturbative QCD predictions introduce 
uncertainties due to missing higher orders. The 
conventional method for its estimation is to con- 
sider the effect of a renormalisation scale vari- 
ation. This method however has at least two 
drawbacks: Since the scale error is positively cor- 
related with a s , this definition produces a bias 
towards small a s values. Secondly there are in- 
dications, that observables calculated from one 
hemisphere only (like the heavy jet mass or -B max ) 
yield less reliable results in the resummation of 
leading logarithms [5]. This should be reflected 
in their theoretical errors. In contrary the scale 
variation yields the smallest uncertainty for the 
heavy jet mass and especially B max . For this rea- 
son a new definition of the theoretical uncertainty 
for the logR prediction was developed. It is based 
on the additional ambiguity which enters the re- 
sumed prediction from the phase space restric- 
tion. These calculations are forced to vanish at 
the kinematic limit y max using the replacement: 
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Usually X = 1 is chosen, although different values 
for this X scale introduce formally only subdom- 
inant contributions [6]. The theoretical error is 
now defined as half of the difference when X is 
varied between 2/3 and 3/2. The same definition 
of the theoretical uncertainty has been adopted 
for the pure NLLA prediction. 

For the 0{a 2 s ) calculation the error from the 
variation around the experimentally optimized 
scales is still used, as in previous publications 

[4 ' 7] ' 

In order to avoid the above mentioned bias, the 
corresponding scale variations are calculated from 
theory distributions with a fixed value of a s . 

3.2. Combination of the results 

For a combination of the a s results from dif- 
ferent observables at the same energy a proper 
treatment of the correlation is mandatory. An 
average value y for correlated measurements yi 



N 



y — Wiyi with: w t 



i=l 



The covariance matrix V has an additive struc- 
ture for each source of uncertainty: 



Ascalc 



V = Y stat _|_ ysys.cxp. _|_ yliad _|_ y 

The statistical component was estimated with 
Monte Carlo simulation, while the correlation of 
systematic errors is modeled by the minimum 
overlap assumption: 



V«=min(A?,A?) 



(2) 



Here Aj denotes the corresponding error of the 
observable i. 

4. Determination of a s from mean values 
with power corrections 

The analytical power ansatz is used to deter- 
mine a s from mean event shapes. This ansatz 
provides an additive term to the perturbative 
0(a 2 3 ) QCD prediction [8]. 



(/> = — / fj-da = </ pert ) + </ pow > 
otot J da 



(3) 



Here an additional non-perturbative parameter 
ao enters for the contributions to the event shape 
below an infrared matching scale /x/. In order to 
measure a s from individual high energy data the 
parameter a has to be determined by a global 
fit of a s and ao to a large set of measurements at 
different energies. 

5. Combination of all DELPHI a s mea- 
surements 

Assuming the validity of the QCD predicted en- 
ergy dependence of a s , all results can be evolved 
to a reference energy, e.g. Mz, and combined to a 
single a s (Mz) measurement. We include also the 
LEP2 results at 133, 161 and 172 GcV from [7]. 
For a s at and around Mz we have reanalyzed the 
distributions from [7] for our five observables and 
combined the results using the same treatment 
of correlation as described in section 3.2. For 
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theory 


a s {M z ) 


exp. 


theo. 


0(aj) 


0.1157 


0.0018 


0.0027 


NLLA 


0.1093 


0.0023 


0.0051 


C(q2) + nLLA 


0.1205 


0.0021 


0.0050 


means 


0.1184 


0.0009 


0.0032 



a (M 7 ) from distributions in 0(a 2 ) 



Table 1 

Results of combining all DELPHI a s measure- 
ments at LEP1 and LEP2. The quoted errors are 
the experimental (stat. and sys.exp.) and theo- 
retical (scale and hadronisation) ones. 



a s from mean values the measurements of events 
with reduced centre-of-mass energy between 44 
and 76GeV [3] have been included as well. In 
combining the a s results one faces again the prob- 
lem of correlation. Although the measurements 
at different energies are obviously statistically in- 
dependent, the systematic and theoretical errors 
are not. Again this part of the covariance matrix 
was modeled assuming minimum overlap. The 
only difference to the method described in sec- 
tion 3.2 is, that the error of ISR and four fermion 
background correlates only those measurements 
which have this source of uncertainty in common. 

The result of these combinations are given in 
Table 1. As an example the Figures 1 and 2 dis- 
play the combination and energy evolution of the 
a s values from the 0(a 2 s ) fit to distributions and 
from mean values with power corrections respec- 
tively. 

As can be seen from Table 1 the total errors 
arc dominated by the theoretical uncertainty. 

6. The running of a s 

In the previous section we have assumed the 
validity of QCD in order to combine the mea- 
surements at different energies. Now we will turn 
around the perspective and use the determination 
of a s at different energies to test the predicted 
scale dependence of the coupling. The logarith- 
mic energy slope of the inverse coupling is (in 
second order) given by: 
-i 

26 + 2bia s - 



89GeV 
91 GeV 
93 GeV 
133 GeV 
161 GeV 
172 GeV 
183 GeV 
189 GeV 
192 GeV 
196 GeV 
200 GeV 
202 GeV 
205 GeV 
207 GeV 



( 0.09) 
( 0.38) 
( 0.09) 
( 0.07) 
( 0.02) 
( 0.01) 
( 0.06) 
( 0.13) 
(0.03) 
( 0.05) 
( 0.02) 
( 0.02) 
(-0.01H 
( 0.03) 



unweighted mean 
weighted funcorr.) mean 
correlated mean 

PDG world average 

j i , \S i i 3 



DELPHI 



dlnQ 



1.27 (QCD) 
0.89 (QCD+gluinos) 



a s (M z ) 

Figure 1. Values for a s at different energies from 
a fit of the 0{a 2 s ) prediction to distributions. Also 
shown are the results of different mean value defi- 
nitions. The numbers in brackets give the weights 
of the individual measurements in the calculation 
of the correlated average. 



Table 2 gives the slopes when fitting the function 
1/(6 log -y/s + c) to the a s values. Figure 2 dis- 
plays the fit result for the a s values from mean 
values with power corrections as an example. The 
results are in good agreement with the QCD ex- 
pectation. The unweighted mean of the different 
measurements yields 1.28 ± 0.12(RMS). In this 
fit the correlation between the measurements is 
taken into account by including the full covari- 
ance matrix into the definition of the \ 2 function. 
The correlation is modeled similarly as described 
in section 3.2. The only differences here are that 
the statistical errors are uncorrelated, and that 
the errors from ISR and four fermion background 
correlate only those data which have these sources 
of uncertainty in common. 

7. The RGI method 

The p functions governs not only the scale de- 
pendence of a s , but directly the energy evolution 
of e.g. event shapes. In the renormalisation group 
invariant (RGI) perturbation theory [1] the en- 
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theory 




dloxEcM ± Stat ± S ^ S 


O(a'i) 




1.27 ± 0.15 ± 0.33 


NLLA 




1.40 ± 0.17 ± 0.44 




fNLLA 


1.32 ± 0.11 ± 0.27 


means 


(pow.corr.) 


1.11 ± 0.09 ± 0.19 



Table 2 

Results for the slope b when fitting the function 
l/(&logy / i + c) to a s values obtained for the dif- 
ferent energies. 



ergy evolution of fully inclusive quantities 1Z, like 
mean values of event shape distributions, can be 
used to measure the [3 function directly: 



dR- 



dQ 



2 " 



fh 



R- 



4.14 (QCD) 

2.76 (QCD+gluinos) 



Here R is 7Z/A, with A the first order coefficient 
when expanding 1Z in powers of a s /n. In [3] 
it was shown, that this approach yields consis- 
tent results for various observables. However, the 
highest accuracy can be achieved for the observ- 
able 1Z = (1 — T) , since here the most low energy 
data is available. Our measurement yields: 



dR- 



dQ 



4.21 ±0.18 



The fit result is displayed in Figure 3. This is in 
good agreement with the QCD expectation and 
allows to exclude the existence of light gluinos in 
the mass range of <30 GeV. 
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Figure 2. Scale dependence of a s from mean val- 
ues with power corrections. The band displays 
the QCD energy evolution of the mean value of 
this measurements. The dotted line indicates the 
fit of the scale dependence. 
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Figure 3. Scale dependence of 1/(1 — T) together 
with theoretical predictions in the RGI frame- 
work and the fit result. 



